Statistical Theory of Energy Transfer to Small and Chaotic Quantum Systems 
Induced by a Slowly- Varying External Field 
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We study nonequilibrium properties of small and chaotic quantum systems, i.e., non-integrable 
systems whose size is small in the sense that the separations of energy levels are non-negligible as 
compared with other relevant energy scales. The energy change AE induced by a slowly-varying 
external field X(t) is evaluated when the range of the variation AA is large so that the linear response 
theory breaks down. A new statistical theory is presented, by which we can predict { AE) , the average 
of AE over a finite energy resolution SE, as a function of X(t) if we are given the density of states 
smeared over SE, the average distance of the anticrossings, and a constant K < 1. 
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Recently, active research has been devoted to small 
and chaotic quantum (SCQ) systems, i.e., non-integrable 
systems whose size is small in the sense that the sepa- 
rations of energy levels are non-negligible as compared 
with other relevant energy scales JlJ]. Such systems are 
frequently encountered in many fields of physics jlj , such 
as excited states of molecules Q and those of quantum 
dots H . SCQ systems have the universal properties that 
level spacings obey a universal distribution, and that the 
energy levels exhibit anticrossings as a function of an ex- 
ternal parameter A Note that these are static prop- 
erties. In contrast, only limited knowledge has been ob- 
tained for dynamical or nonequilibrium properties of SCQ 
systems A major difficulty is that solving the time- 
dependent Schrodinger equation (TDSE) becomes impos- 
sible when / > 4, where / is the degrees of freedom of 
the system, for the most common case where the "coordi- 
nate" of each degree of freedom takes continuous values. 
This is because the computational time generally grows 
exponentially with /. If / were huge enough, on the 
other hand, the nonequilibrium thermodynamics would 
be applicable, by which one could predict nonequilibrium 
properties much more easily than by solving the TDSE. 
For SCQ systems, however, the nonequilibrium thermo- 
dynamics is not applicable because its basic assumptions, 
such as the local equilibrium, are not satisfied. The pur- 
pose of the present Letter is to propose a new theory, 
called the random-probability netweork (RPN) method, 
by which nonequilibrium properties of SCQ systems can 
be predicted easily, without solving the TDSE ||. 

We consider an SCQ system subject to an external field 
A(t), which represents any external disturbance, such as 
a magnetic or electric field, and the position coordinate 
of a moving "wall" (a portion of another system which 
couples to the SCQ system). The energy change AE 
induced by the variation of A(t) |(| is evaluated under 
the conditions that (i) the range of the variation AA is 
large so that AA cannot be treated as a small perturbation 
and the linear response theory breaks down, and that (ii) 



the typical time scale T of the variation of A is long (so 
that r of Eq. @ < 1). Without solving the TDSE, we 
can easily predict (AE 1 ), the average of AE over a finite 
energy resolution 8E, by the RPN method if we are given 
a constant K (< 1) and two "fundamental functions" 
which contain only small information of the system. 

To be concrete, we explain the RPN method using a 
simple model. However, since the basic idea is based on 
the universal properties mentioned above, we expect that 
the RPN method is generally applicable to SCQ systems. 
Consider the time evolution of coupled rotators (angular 
momenta L and M), whose Hamiltonian is @] 



H(t) = A(t) ( 



n 1 (t)+n 2 , (i) 



in an appropriate unit system, in which the energy, A 
and h are dimensionless The external field X(t) is 
taken as X(t) = 1 - AAjl - (£/T) 2 ] 4 for -T < t < T and 
A(t) = 1 for \t\ > T S. Since we are interested in the 
case where AA is large, we here take AA = 0.5. For this 
X(t), 7i takes the same form at both ends, t — ±T, and 
AE can be finite only when transitions occur between 
different levels. However, it should be stressed that the 
RPN method is also applicable to the case of 7i(T) ^ 
"H(—T), for which (AE) consists of both the shifts of 
energy levels and the transitions between them. 

As in the case of the statistical mechanics, we must first 
classify quantum states according to constants of motion, 
to obtain a set of subspaces, each of which has no con- 
stant of motion. It is in each subspace that the system has 
chaotic natures and the RPN method is applicable. We 
here consider the subspace of J z /h — even, I = m = 15, 
E = +. Here, J z = L z + M z , h 2 l(l + 1) and h 2 m(m + 1) 
are the quantized values of L 2 and M 2 , and E denotes 
the parity under the exchange of L z and M z . We take 
h = 0.322748612, so that the corresponding classical sys- 
tem is chaotic 0. The eigenvalues E n (X) exhibit anti- 
crossings as functions of A as shown in Fig. [l| In the inset 
of Fig. ^|, we plot the density of states p$e(E; A) which is 
smeared over a finite energy range SE. Here and after, 
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we consider the E < part only, because the system is 
symmetric about E = 0. On an average, pse(E; A) is an 
increasing function of E for E < (whereas it is decreas- 
ing for E > 0), and the average energy spacing (= 1/pse) 
decreases with E. We have confirmed that the normal- 
ized spacing, s„(A) = pse(E; X)(E n+1 (X)-E n (X)), obeys 
the Wigner distribution, P(s) ~ (tts/2) exp(— 7rs 2 /4) [Q. 
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FIG. 1. The energy levels as a function of A. 
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FIG. 2. The integrated energy change Ase(E) for T — 
100, 150, 200, obtained by solving the TDSE (solid lines) 
and by the RPN method (dashed lines). Inset: The smeared 
density of states pse(E; A) for A = 1.0, 8E — 1.0. 

Before presenting the RPN method, we first explore 
dynamical properties of the system by numerically solv- 
ing the TDSE, lU§i\^{t)) = H(t)\ip(t)). Here, we use the 
split operator method || to obtain the time evolution 
operator U(t 2 ,h) defined by \ip(t 2 )) = U{h,t x )\^{h))- 
Noting that li{t) is constant for \t\ > T, we take one 
of its eigenstates \ip n ), where Ti.(\t\ > T)\ijj n ) = E n \ip n ), 
as the initial state. The state at t > —T is denoted by 



\ifj[ n ](t)} = U(t, —T)\ip n ). We are interested in the energy 
change AE n between the initial (t = — T) and the final 
(t = T) states §; AE n = ^ [n] (T)\n(T)\^T)) - E n . 
This is plotted against E n for T = 100 in Fig. |[ It is seen 
that AE n scatters rather randomly as a function of E n . 
It is therefore impossible to predict AE n of individual 
levels, without solving the TDSE. 
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FIG. 3. The energy change for T — 100, as a function of 
the initial energy E n . 

To see what predictions are possible and meaningful, 
let us consider typical experiments on SCQ systems. Sup- 
pose that one tries to prepare the initial state in \ip n ). 
However, except for the ground state, it is difficult to pre- 
pare the initial state of SCQ systems accurately. Hence, 
one can actually prepare the initial state in a finite energy 
range E n — 8E/2 < E < E n + 5E/2, where SE denotes 
the energy resolution of the apparatus. Therefore, the 
initial state should be an incoherent mixture or a coher- 
ent superposition of |^ m )'s of \E m — E n \ < SE/2. Let 

AEgg(E) be the measured value of the energy change 
in the jth one of such experiments, where E = E n . 
This quantity fluctuates from experiment to experiment 
for two reasons: One is the random fluctuation of AE n 
as a function of E n mentioned above; this fluctuation 
appears in AE^l (E) because the intial state is prepared 
randomly in the energy range \E m — E\ < 5E/2. The 
other reason is the quantum-mechanical fluctuation; even 
if 5E — > 0, only the average of AEg^,(E) over many ex- 
periments would agree with AE n . Namely, even for the 
same initial state \ipn), the measured value of each exper- 
iment scatters about AE ni with the standard deviation, 
D n ee [(ip[ n ](T)\(H(T) - E n ) 2 \ip [n ](T)) ~ (AEn) 2 } 1 / 2 . By 
numerically computing D n , we find that D n ~ AE n . 
Because of these large fluctuations, a single experiment 
is insufficient for comparison of the experimental result 
with a theoretical one, for an SCQ system. We must 
therefore perform many experiments, and we are most 
interested in the average AEse(E) of their results. We 
will argue later, and have confirmed numerically, that 
the quantum coherence of the initial state is irrelevant 
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to this quantity, under our assumption that AA is large 
so that the system experiences many anticrossings during 
the variation of X. Therefore, AEg E (E) equals the aver- 
age of AE m 's over m's which satisfies \E m — E\ < SE/2] 



AE SE {E) 



E+SE/2 



E-&E/2 



de ^E n 



b& E (e - E m ) 



E+SE/2 
E-SE/2 



(2) 



de^tpSE^ - E m ) 



where 4>se(&) is the normalized gaussian function with 
width SE. To find statistical tendencies of AEse{E), we 
investigate the integrated energy change defined by 



Ase(E) 



de-^2AE n cj) SE (e- E„ 



(3) 



We can obtain AEse(E) from Ase(E) as 
A 5E {E + SE/2) - A SE {E - SE/2) 



Pse{E)5E 



AE SE {E). (4) 



We plot As E (E) for SE = 1.0, as a function of E, by the 
solid lines in Fig. ||. We find that A$ E (E) is a smooth 
and, on an average, increasing function of E. It strongly 
suggests that the prediction about As E (E) may be pos- 
sible by a simple theory, without solving the TDSE. 
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FIG. 4. A pair of the RPN diagrams (a) and (b). 

As such a theory we now present the RPN method, 
by which we can calculate A$ E {E) (hence (AE)) as a 
function of A(i), if we are given a constant K and two 
"fundamental functions," which are ps E (E;X) and the 
average distance X S (E; A) (along the A axis of Fig. |l]) of 



the anticrossings [|10|. Note that these functions contain 
only small information on the system because they rep- 
resent the distribution of energy levels only, and further 
loss of information occurs by averaging it over a finite en- 
ergy SE. The RPN method enables us to predict As E (E) 
from such small and reduced information. 

The RPN method is applicable under the conditions 
that (i) AA is large so that each level encounters many 
anticrossings from t = —T to T, and that (ii) T is long 
so that most transitions between different levels occur 
only at anticrossings, and the transition at each anti- 
crossing occurs only between adjacent levels. (This is 
satisfied if r of Eq. (||) <C 1.) Both conditions are sat- 
isfied for A A = 0.5 and T > 100 for the present model. 
Under these conditions, anticrsossings may be regarded 
as "bridges," only through which transitions occur, and 
the RPN method is formulated as follows. 
Step 1: Associate the level diagram of the original sys- 
tem, Fig. [l], with many pairs of hypothetical diagrams, 
which we call the RPN diagrams. Each pair consists of 
two diagrams (a) and (b), as shown in Fig. ||. Diagram 
(a) is a set of hypothetical energy levels • • • , E%-, 

where the number N of the levels are taken equal to 
or larger than the number of relevant energy levels. In 
each small interval [A, A + dA), the hypothetical levels are 
generated randomly according to ps E (E;X) of the origi- 
nal system and the Wigner distribution of the spacings. 
Here, the small interval dX is taken in such a way that 
dX/X s (E*(X); A) <C 1. Diagram (b) is a set of straight 
lines (labeled by n — 1,2, •••) with "bridges" (labeled 
by B\, B 2 , • • •)• The bridges are randomly generated in 
such a way that the probability of finding one bridge 
in each small interval [A, A + dX) on line n is given by 
dX/X s (E*(X); A), where E*(X) is given by diagram (a). 
We also impose the condition that any two bridges are 
not on a single line in the same interval. This constric- 
tion is consistent with the assumption dX/X s <C 1, which 
states that the probability of finding such two bridges in 
the original system is negligible. 

Step 2: Attach an NxN matrix w to each bridge, where 
w(B) at a bridge B between lines n and n + 1 is given 
by w n n+ i = w n+ i „ = r, w n . n = w„+i,„+i = 1 - r, and 



Si j for i ^ n,n+ 1. Here, 



r((E* 



E* l+1 )/2;X) represents the transition probability at the 
bridge, where 

r(E; A) = exp [-2nK 2 X s (E; X)/h\X\p SE (E; A)] . (5) 

Here, K is a constant (< 1), which is considered as a 
fundamental parameter of the RPN method. This form 
of r(E; A) is deduced from the Landau-Zener formula for 
two levels; r = exp(— 2ttV 2 /h\X\\Fi — F 2 \). Here, by re- 
placing the difference \Fi — F 2 \ between the slopes of the 
levels with K 2 / pse(E; X)X S , and the energy gap V at 
the anticrossing with K^/ pg E (E; A), where K 2 — 1 and 
K 3 < 1, we arrive at Eq. (|) with K = K z j^/W 2 < 1. 
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Step 3: Construct an N x N matrix W as follows. Draw 
the vertical line at A = A(— T) in diagram (b). Asso- 
ciate the vertical line with the matrix W, and take it 
as the unit matrix. Move the vertical line as a func- 
tion of t, according to A = X(t). When the vertical 
line encounters a bridge B, multiply W by w(B) of the 
bridge: W — * w(B)W. If the line encounters two or 
more bridges simultaneously, multiply W by w's of the 
bridges, where, the order of w's are arbitrary because 
such w's of the bridges on the same vertical line com- 
mute with each other. In this way, W at the final time 
(t = T) becomes the product of w(B)'s along the "path" 
of A(t). We interpret its element W mn as the probabil- 
ity of the transition from line n at t = —T to line m at 
t = T, of the random probability network. 
Step 4: Evaluate the average change of the hypotheti- 
cai energy by AE* = £ m (££(A(T)) - E*(\(-T))W mn , 
which is the counterpart of AE n of the original system. 
Calculate the integrated energy change of the RPN by 

r E 

A* SE (E) = de-J2 ^KMe - E*). (6) 

J-oo n 

Step 5: Repeat steps 1-4 many times, to obtain many 
Ag E (E)'s, which correspond to different pairs of the RPN 
diagrams. Calculate the average of Ag E (E)'s. This aver- 
age is the result for Ase(E) by the RPN method. 

The total computational time is tremendously reduced 
as compared with solving the TDSE; e.g., from a few days 
to a few seconds for the model of Eq. ([!]) . The reduction 
becomes much more drastic as / is increased. 

To check the validity of the RPN method, we compare 
Ase(E) of the RPN method with that obtained by solv- 
ing the TDSE. In the RPN method we use p SE (E;X) 
and A S (-E;A) that are obtained by solving the time- 
independent Schrodinger equation. The constant K is 
determined so as to optimize the agreement between the 
results of the RPN method and the TDSE. (Hence, K is 
the only free parameter in the comparison.) This gives 
K = 0.177, in consistent with K < 1. The dashed lines in 
Fig. plot Ase (E) obtained by the RPN method for vari- 
ous values of T. On an average, they agree with the solid 
lines obtained by solving the TDSE, although the fine 
structures, peaks and dips, of the solid lines are absent 
in the dashed lines. It is expected that such peaks and 
dips would be shallower for larger systems. Considering 
also the tremendous reduction of the computaional time, 
the overall agreement seems satisfactory. The reason why 
the RPN method gives good results is partly that we have 
introduced the finite energy resolution 5E, which results 
in the average over many initial states. This smears out 
strong dependence on the initial state. Another reason is 
that we have assumed a large value of AA. As a result, 
the wavefunction experiences many anticrossings under 
the variation of A(t), and effects of the phase coherence 
on (AE) are almost destroyed on an average. Detailed 



discussions on this point will be described elsewhere. 

As we have done above, we may obtain the fundamen- 
tal functions by solving the time- independent Schrodinger 
equation, which is much easier than solving the TDSE. 
For larger systems, or for systems whose Hamiltonian 
is unknown, the fundamental functions and K may be 
obtained from some experiments (e.g., on the specific 
heat) and/or experiences (in similar systems), as in the 
case of estimating thermodynamical functions from ex- 
periments. To develop a systematic method to perform 
this program is a subject of future study. 

We finally note that Ase(E) is, on an average, an in- 
creasing function of E. This means that if X(t) is varied 
repeatedly, then, as a result of accumurated transitions, 
the energy of the system tends toward E — 0, where 
Pse(E) is highest. We thus find that the system tends to 
"climb" the curve of the smeared density of state. There- 
fore, an SCQ system tends to absorb the energy if its 
Pse(E) is an increasing function of E. 
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